Negative differential conductance in quantum dots in theory and experiment 



M. C. Rogge 1 , F. Cavaliere 2 ' 3 , M. Sassetti 2 , R. J. Haug 1 , and B. Kramer 3 

1 Institut fur Festkorperphysik, Universitdt Hannover, Appelstrafie 2, 30167 Hannover, Germany 
2 /. Institut fur Theoretische Physik, Universitdt Hamburg, Jungiusstrafie 9, 20355 Hamburg, Germany 
3 INFM-LAMIA, Dipartimento di Fisica Universita di Genova, Via Dodecaneso 33, 16 146 Genova, Italy 

(Dated: February 2, 2008) 

Experimental results for sequential transport through a lateral quantum dot in the regime of spin 
blockade induced by spin dependent tunneling are compared with theoretical results obtained by 
solving a master equation for independent electrons. Orbital and spin effects in electron tunnel- 
ing in the presence of a perpendicular magnetic field are identified and discussed in terms of the 
Fock-Darwin spectrum with spin. In the nonlinear regime, a regular pattern of negative differential 
conductances is observed. Electrical asymmetries in tunnel rates and capacitances must be intro- 
duced in order to account for the experimental findings. Fast relaxation of the excited states in 
the quantum dot have to be assumed, in order to explain the absence of certain structures in the 
transport spectra. 

PACS numbers: 73.63.Kv, 71.10.Pm, 72.25.-b 



Since the discovery of the Coulomb blockade effect in 
semiconductors 1 electronic transport properties of quan- 
tum dots have been the subject of continuous interest. 
The interplay of Coulomb interaction between quantized 
charges, the spectra of orbital states and the electron 
spin causes a rich variety of effects in the transport spec- 
tra of lateral quantum dots 2 . In an external magnetic 
field, signatures have been found for transitions between 
different angular momentum states of the single elec- 
tron Fock-Darwin spectrum 3 ^, modified by interaction 
effects 5 . Effects of the electron spin have been observed, 
including several transport blocking mechanisms. One of 
these spin blockades was introduced by Ciorga et al££. 
Due to spin polarized leads Coulomb peaks show a spin 
dependent amplitude. This was observed in split gate 
quantum dots and recently in devices made with local 
anodic oxidation 8 . In nonlinear measurements negative 
differential conductance was observed and explained in 
terms of spin dependent tunneling 9 -. 

In the present paper, we report results of spin block- 
ade measurements of the linear and nonlinear sequential 
transport of a 2D quantum dot in a perpendicular mag- 
netic field in comparison with a theoretical model. In 
the considered parameter region, we confirm results con- 
cerning angular momentum and spin blockade effects ob- 
tained earlier 6,7,9 . In particular, we observe the charac- 
teristic zig-zag patterns of the linear conductance peaks 
when varying the magnetic field. The positions and the 
strengths of the linear conductance peaks show a reg- 
ular alternating behavior as a function of gate voltage 
and magnetic field. In the nonlinear transport regime a 
complicated substructure involving certain transitions to 
excited states is revealed. A striking bi-modal pattern 
of positive and negative differential conductance (PDC 
and NDC) is observed below the (N + i) —> N transition 
line which is the extension of the linear zig-zag pattern. 
This appears to be associated with a strong depletion 
of the conductance peaks corresponding to transitions 
N — > (N + 1). To study the conditions for obtaining 
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FIG. 1: (a) AFM picture of the gate structures used to define 
electrostatically a lateral quantum dot in the inversion layer 
of a AlGaAs/GaAs heterostructure. (b) electrostatic model 
for the device in a perpendicular magnetic field. The dot 
with edge and core is coupled via tunnel resistances R and 
capacitances C to source S and drain D that feature spin split 
edge channels. The in-plane gate G is coupled capacitively. 



the NDC in the nonlinear transport induced by the spin 
blockade we compare the experimental findings with the- 
oretical results for the nonlinear current- voltage charac- 
teristics that were obtained by solving the rate equation 
for the Fock-Darwin spectrum with suitably defined spin 
dependent tunnel rates. In addition to the linear con- 
ductance we investigated in detail the conditions for the 
spin induced NDC to appear in the nonlinear differential 
conductance. It is shown that assuming spin dependent 
tunneling only is not sufficient to understand the data. 
Our main result is that it is necessary to assume a rather 
complex interplay of asymmetries, and the relaxation of 
certain excited states to be much faster than the transi- 
tion rates between leads and core states of the quantum 
dot, in order to account for the regularity in the exper- 
imental conductance traces in a large parameter region. 

The quantum dot has been fabricated in an inversion 
layer embedded in a GaAs/AlGaAs heterostructure. The 
2D electron system (2DES) is located 34 nm below the 
surface. The sheet density is n = 4.3 • 10 15 m 2 . A combi- 
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nation of Local Anodic Oxidation (LAO) with an Atomic 
Force Microscope (AFM) and electron beam lithography 
(ebeam) is used to define the quantum dot structure vis- 
ible in Fig. QJl. Details of the sample processing can be 
found in [10,11]. Quantum dot A is coupled to source 
(S) and drain (D) leads via resistances Rs and Rd and 
capacitances Cs and Cd (Fig. Wp)- An in-plane gate G 
is coupled via Cg and can be biased to tune the tun- 
neling rates and electron numbers of the dot. Although 
the dot B is coupled capacitively to A the measurements 
discussed here reflect the dynamics of dot A only. They 
are performed in a regime with dot B only connected 
to the source lead. Thus transport can occur only via 
dot A. The influence of the second dot is discussed in 
[8] . With dot A in the Coulomb blockade regime the dif- 
ferential conductance is recorded using standard lock in 
technique at a temperature of 50 mK. 

Figures show typical results for the differential 
conductance in the plane of Vq and B with V = 0. The 
sawtooth-like traces that are separated by the charging 
energy show several characteristic features. Firstly, there 
is a clear bimodal behavior of the peak position. An os- 
cillation with a huge increase follows one with a small 
increase. The same is true for the intensity of the con- 
ductance due to spin dependent transport. Secondly, the 
intensity of all traces with positive slopes is very weak, al- 
most not visible in this plot. Thirdly, traces correspond- 
ing to successive Coulomb peaks are shifted relative to 
each other by half a period. Large intensities in one trace 
correspond to small and weak features in the neighboring 
traces. In the nonlinear conductance (Fig. Efc,d), such 
features are also present. Additionally, one notes that 
while the large and strong features are associated with 
positive differential conductance in the nonlinear region, 
the smaller and weaker features are continued into the 
nonlinear regime with NDC. An important observation 
is that only traces with a negative slope as a function of 
the magnetic field are observed to yield NDC. 

In order to understand these features, we consider a 
model for N uncorrelated electrons in a 2D harmoni- 
cally confined quantum dot. The ground state energy 
for N electrons in this system is Hn — EcN 2 /2 — 
eNV G C G /C^ - eNVCs/C^ + £f =1 ^(£). The first 
term is the classical charge addition energy (Ec charg- 
ing energy) which accounts for the Coulomb repulsion. 
Cs = Cs + Cd + Cg ~ Cs + Cd is the total capaci- 
tance of the dot. The lever arm a = Cq/C^ accounts 
for the energy change induced by Vq- The last term is 
the kinetic energy contributed by occupying the ener- 
getically lowest one-electron states of the quantum dot 
with N electrons. The chemical potential for the transi- 
tion N — 1 <-> N is defined as Hn(B) '•= Hn — H^-i = 
E C (N - 1/2) + E N (B) - eV G C G /C^ - eVC s /C^ For 
the nonlinear regime also the lowest lying excited states 
with energies = Hn + en+i — £n must be considered. 

For harmonic confinement, £j(B) = Ej(B) + SjE z (B) 
with E z = g* hbB. The first term represents the j-th 
level of the Fock-Darwin (FD) spectrum, while the second 
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FIG. 2: (a) Linear conductance peaks of the quantum dot 
in the plane of gate voltage Vg and magnetic field B (dark: 
high positive differential conductance, bright: low conduc- 
tance); (b) enlarged inset of (a). The Coulomb peaks develop 
a bimodal behavior in position and amplitude, (c) spectrum 
of the nonlinear conductance in the same parameter region 
at bias V = 0.5 mV (red: NDC); (d) enlarged inset of (c). A 
regular pattern of PDC and NDC is observed. 



term describes the Zeeman splitting of the FD spectrum 
with g* the electron ^-factor, the Bohr magneton 
and Sj = ±1 the z component of the spin of the j-th FD 
level (units H/2). Here, we assume S\(B) < 62(B) < — 
For small B, the ground state is a spin singlet or triplet 
depending on whether N is even or odd, respectively. In- 
creasing the magnetic field, many level-crossings occur 
in the FD spectrum. This is shown in Fig. where the 
FD spectrum and the energy level E^(B) are shown for 
dot filling factor (u) 2 < v < 4. For very large magnetic 
field, uo c = eB /m* ^> c^o (m* electron effective mass, c^o 
characteristic frequency of harmonic confinement poten- 
tial), the Fock-Darwin levels merge into almost degener- 
ate, Zeeman split Landau bands. 

The dot is connected to external noninteracting leads 
by means of tunneling junctions with resistances Rs,b 
and capacitances Cs,d (Fig- \3p)- The source (S) ter- 
minal is kept at a potential V, while the drain (D) is 
grounded. At low temperature, for V > electrons flow 
essentially from the drain to the source lead. We consider 
the tunneling in the sequential regime and introduce the 
tunneling rates ex Jij^o with To in dimension of a 

tunneling rate and jij dimensionless parameters which 
depend on the states involved in the tunneling process. 
In particular, we have found to be useful to introduce the 
following parameters: js,D tunneling via source or drain; 
7 e?c tunneling via edge or core ground state transitions; 
7* c tunneling via edge or core excited states transitions; 
7l ? l tunneling with spin down or spin up. For example, 
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for tunneling with spin up via the source to an excited 
state in the core we have 7&,/ = 7f • 7s • 7* (throughout 
this work, j e = 1, 7* = 1 and 7^ = 1). 

We have calculated the current and the differential con- 
ductance by solving a master equation. First, we address 
the linear transport regime where only ground states TV, 
TV + 1 are involved. The transition TV <-> TV + 1 occurs if 
the resonance condition /ijv+i(5) ~ is fulfilled. There- 
fore, the peak position in the plane {B,Vq) reflects the 
position as a function of B of the TV + 1-th energy level. 
By comparing the linear conductance traces in Fig. [2Jl 
and the number of kinks in the trace of the chemical po- 
tential in Fig. one can conclude that the number of 
electrons is in the range TV = 53, . . . , 57. In the 2 < v < 4 
regime, the dot states with an increasing (decreasing) 
Ej (B) have a low (high) angular momentum and are lo- 
cated in the core (edge) of the dot (Fig. QJ)). This leads 
to a weaker (stronger) overlapping of the electron wave- 
function with the leads states. To take into account this 
fact, tunneling rates involving a low angular momentum 
state are characterized by the factor < j c < 1. More- 
over, because of the spatial separation of spin-resolved 
edge states in the leads, tunneling processes involving a 
state with spin Sj = +1 are depressed. We assume that 
the corresponding rates are renormalized by < 7^ < 1. 
With these assumptions, the bi-modal behavior can be 
explained. As shown in Fig. 03,c where the linear con- 
ductance is plotted, lines with positive slopes are strongly 
suppressed because of the tunneling involving states in- 
side the core of the dot. The lines with a negative slope 
exhibit a pronounced bi-modal behavior in the intensity, 
which reflects the alternating spin of the TV+l-th electron 
in the FD spectrum, Fig. OK 

At first glance, it might seem that the above assump- 
tions concerning the tunneling rates are quite exhaus- 
tive, and should be sufficient to understand also the non- 
linear behavior at least qualitatively. However, this is 
not the case, as we have verified by solving the master 
equation including the states relevant to the nonlinear 
regime. These are the ground states for TV and TV + 1 
electrons as in linear transport and their lowest-lying ex- 
cited counterparts TV* and TV + 1*. For V > 0, due to 
the difference in the chemical potentials of source and 
drain, the TV <-> TV + 1 conductance trace splits, giving 
rise to "stripes" in the (£>, Vq) -plane (Fig. EK). In addi- 
tion, new traces due to the transitions TV <-> TV + 1* and 
TV* <-> TV + 1 appear. The tunneling rates involving an 
excited state of the dot with low angular momentum are 
renormalized by the factor < 7* < 1. 

To obtain NDC in the upper part of the stripes 
(Fig. Eli, transitions TV + 1 -> TV and TV+ 1 -> TV*) a tun- 
neling asymmetry 75/71) > 1 must be assumed, so that 
the state TV* is dynamically trapped in the do1iA2iAiAi and 
the current decreases. However, this yields as a function 
of B traces of conductance peaks with NDC which have 
both, positive as well as negative slopes (Fig. 0t>). In 
contrast to the experimental findings, the NDC features 
with positive slopes are even more pronounced than the 
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FIG. 3: (a) Fock-Darwin energy spectrum of a quantum dot 
as a function of the magnetic field B in the 2 < v < 4 regime: 
states with spin —1 (+1) correspond to solid (dashed) lines; 
bold: E 57(B). (b) Theoretically evaluated linear conductance 
traces as a function of Vb and B for the transitions TV = 
52 -> 53,. . . ,N = 56 -> 57 with 7t 0.5 and 7 C 0.2 (see 
text), (c) Enlarged region of (b); parameters: hu>o = 1.3 meV, 
g * 0.44, T = 70 111K. 



ones with negative slope. The tunneling suppression due 
to dot states with low angular momentum enhances the 
dynamical trapping thus favoring NDC. Moreover, in ex- 
periment the transition lines corresponding to transitions 
TV — > TV + 1 and TV — > TV + 1* are clearly suppressed. 

For a description consistent with the experiment, pa- 
rameters have to be chosen such that the NDC traces 
with positive slopes vanish. In addition, the intensity of 
the lower transition lines must vanish. The former can be 
achieved assuming 7* < j C: and by assuming a strong re- 
laxation of the excited states towards the corresponding 
ground states. We have introduced a phenomenological 
relaxation rate tTq (r dimensionless constant) into the 
master equation and consider r ^> 7*. Then, the relax- 
ation time is much faster than the dwell time of the low 
angular momentum excited states, leading to a substan- 
tial reduction of the NDC. The differential conductance 
is shown in Fig.^fc. Due to the strong relaxation the con- 
ductance exhibits spin-driven bi-modal negative-positive 
differential conductances in traces with negative slopes, 
and no NDC in traces with positive slopes. 

In order to obtain vanishing conductance peaks in the 
lower stripes, it is necessary to consider asymmetric ca- 
pacitances. In Fig. EJi the differential conductance for the 
same set of parameters as in Fig. is shown, with an 
asymmetry Cs/Cd = 1/9- Qualitative agreement with 
the experimental findings is obtained. 

We have attempted to reproduce theoretically the ex- 
perimental findings by using many different, physically 
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FIG. 4: (a) Scheme of the nonlinear conductance traces in 
the (B,Vg) plane in the voltage region of the experiment; 
(b) nonlinear conductance for asymmetric tunnel barriers 
Js/jd — 50 and 7* — j c (red: NDC); (c) nonlinear con- 
ductance as in (b) with strong relaxation of excited states 
r = 1.1 and 7* = 0.257 c ; (d) nonlinear conductance as in (c) 
and asymmetric capacitances Cs/Cd = 1/9. 



reasonable, combinations of parameters. However, only 
with the set of parameters discussed here we were success- 
ful. Apparently, the patterns of the nonlinear differential 
conductance traces with the properties shown in Fig. 
depend very sensitively on the properties of the quantum 
dot and its environment, in addition to the properties of 
the tunnel contacts. This has been confirmed experimen- 
tally by measuring the nonlinear differential conductance 
of a different 2D quantum dot. In this case, NDC could 
not be achieved, although the above discussed bi-modal 



features of the linear conductance traces were detected. 

In conclusion, we have compared measurements of the 
nonlinear conductance spectra of 2D quantum dots in 
a strong perpendicular magnetic field with theoretical 
current- voltage characteristics obtained by solving a mas- 
ter equation using the Fock-Darwin model together with 
a set of parameters that characterize the coupling be- 
tween the quantum dot and the leads, and the properties 
of the electron states. We have found that in order to 
reproduce the most prominent features observed exper- 
imentally, a rather complex set of model parameters is 
needed. We have found evidence that it is not sufficient 
to assume spin dependent tunneling. Especially for ex- 
plaining the regular bi-modal patterns of NDC in the 
nonlinear conductance it appears to be necessary to as- 
sume fast relaxation of certain excited quantum dot states 
in addition to asymmetry of the tunnel rates. Further- 
more, for reproducing the experimentally observed strik- 
ing suppression of the traces corresponding to transitions 
N — > (N + 1) and N — > (TV + 1)*, asymmetry in source 
and drain capacitances are necessary. 

In contrast to the Coulomb blockade effect, which is 
universal and independent of the details of the model 
in the region of sequential tunneling, the spin blockade 
effect appears not to be universal but very delicately de- 
pending on details of the experimental setup. In addition, 
the experimental data shown in Fig. exhibit features, 
like the behavior of the slopes of the conductance traces 
with magnetic field, that apparently cannot be under- 
stood within the Fock-Darwin model including the spin. 
For understanding these features, it appears necessary to 
take into account electron correlations beyond the con- 
stant interaction model. 
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